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Abstract 

The optical absorption of a many (continuum) polaron gas is derived in the 
framework of a variational approach at zero temperature and weak or inter- 
mediate electron-phonon coupling strength. We derive a compact formula 
for the optical conductivity of the many-polaron system taking into account 
many-body effects in the electron or hole system. Within the method pre- 
sented here, these effects are contained completely in the dynamical structure 
factor of the electron or hole system. This allows to build on well-established 
studies of the interacting electron gas. Based on this approach a novel feature 
in the absorption spectrum of the many-polaron gas, related to the emission 
of a plasmon together with a phonon, is identified. As an application and 
illustration of the technique, we compare the theoretical many-polaron op- 
tical absorption spectrum as derived in the present work with the 'd-band' 
absorption feature in Nd2Cu02. Similarities are shown between the theoret- 
ically and the experimentally derived first frequency moment of the optical 
absorption of a family of differently doped Nd2-3:Cea,Cu04_jy materials. 
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I. INTRODUCTION 



As is the case for polar semiconductors and ionic crystals insight into the nature 
of polarons in high-temperature superconductors can be gained by studying the optical 
properties of these materials. The goal of the present paper is to present a theory of the 
optical conductivity of a system of continuum polarons at any density, including many- 
body effects between the constituent charge carriers, for small and intermediate values of 
the electron-phonon coupling constant and for zero temperature. The method we develop 
here is based on the variational method introduced by Lemmens, Devreese, Brosens (LDB) 
for the ground state energy of the many-polaron gas. The advantage of this approach 
over other theories of many-polaron optical absorption is that it allows to include the 
many-body effects in the system of the constituent charge carriers of the polaron gas on the 
level of the dynamical structure factor of the underlying electron (or hole) system. Thus it 
is possible to select the level of approximation used in the treatment of the many-polaron 
gas by choosing the appropriate expression of the dynamical structure factor for the electron 
(or hole) system. 

Recently the infrared spectrum of cuprates has been the subject of intensive in- 
vestigations PH13|, especially in the case of the neodymium-cerium cuprate family 
Nd2_a;Ce2.Cu04_y (NCCO) PJ5|-pr5|. Several optical absorption features in the infrared 
cuprate spectrum have been tentatively associated with large polarons [0,0 or with a 
mixture of large and small (bi)polarons [|l^. These comparisons with polaron theory were 
derived using a single-polaron picture, so that the density (doping) dependence of the optical 
absorption spectra could not be studied in detail. The many-body theory of the A^-polaron 
spectrum, presented here, allows to study the density (doping) dependence of optical absorp- 
tion spectra. As a first application of the many-polaron optical absorption theory introduced 
here, a preliminary comparison is presented between the theoretical many-polaron optical 
absorption derived in the current work and the mid-infrared spectrum of the neodymium- 
cerium cuprates recently determined experimentally by Calvani and co-workers 



2 



II. OPTICAL ABSORPTION IN THE MANY-POLARON SYSTEM 



A. The LDB variational wave function for a many-polaron system 
The Hamiltonian of a system of N interacting continuum polarons is given by: 

j=l^"''b k k j=l ^^oo j=l£(^j)=l 1^* ■'il 

(1) 

where r^, pj represent the position and momentum of the constituent electrons (or holes) 
with band mass mi,; a^,ak denote the creation and annihilation operators for longitudinal 
optical (LO) phonons with wave vector k and frequency ci^lo; H; describes the amplitude 
of the interaction between the electrons and the phonons; and e is the elementary electron 
charge. The ground state energy of this many-polaron Hamiltonian has been studied be- 
fore by LDB 0, for weak and intermediate strengths of the electron-phonon coupling, by 
introducing a variational wave function: 

I^^LDb) =t/|0)|</^el), (2) 

where \ipei) represents the ground-state many-body wave function for the electron (or hole) 
system and |0) is the phonon vacuum, and ?7 is a many-body unitary operator which deter- 
mines a canonical transformation for a fermion gas interacting with a boson field: 

U = exp |eE (/k^ke*''-^^ - /k^e-^-^^) | . (3) 

In the limit of one fermion, U reduces to a canonical transformation inspired by Tomonaga 
p!6| and applied later by several workers, after Lee, Low and Pines [0, but always for 
one particle-theories. In LDB [Q, this canonical transformation was extended and used to 
establish a many-fermion theory. The /k were determined variationally P| resulting in 



A = W 

hULO + 



2mbS{k) 



for a system with total momentum P = J2j Pj = 0. In this expression, S'(k) represents the 
static structure factor of the constituent interacting many electron or hole system : 

/ N N \ 

\i=ij'=i / 

The angular brackets (...) represent the expectation value with respect to the ground state. 
It may be emphasized that (|^), although it appears like a straightforward generalization 
of the one-particle transformation in [^], represents -especially in its implementation- a 
nontrivial extension of a one-particle approximation to a many-body system. As noted in 
the introduction, the main advantage of the LDB many-polaron variational approach lies in 
the fact that the many-body effects in the system of charge carriers (electrons or holes) are 
completely contained in the structure factor of the electron (or hole) gas. This advantage 
will be carried through into the calculation of the optical properties of the interacting gas 
of continuum polarons which is the subject of the current paper. 

B. Kubo formula for the optical conductivity of the many-polaron gas 

The many-polaron optical conductivity is the response of the current-density, in the 
system described by the Hamiltonian (|l]), to an applied electric field (along the x-axis) with 
frequency uj. This applied electric field introduces a perturbation term in the Hamiltonian 
(0), which couples the vector potential of the incident electromagnetic field to the current- 
density. As is well known, within linear response theory, the optical conductivity can be 
expressed through the Kubo formula as a current-current correlation function ||18| : 



In this expression, V is the volume of the system, and Jx is the x-component of the current 
operator J, which is related to the momentum operators of the charge carriers: 

J = -Ep. = -P, (7) 
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with q the charge of the charge carriers (+e for holes, — e for electrons) and P the total 
momentum operator of the charge carriers. The real part of the optical conductivity at 
temperature zero, which is proportional to the optical absorption coefficient, can be written 
as a function of the total momentum operator of the charge carriers as follows : 



Re[aH] = Rejy^ e^^* ([P.(t), P.(0)]) rftj . 



(8) 



Integrating by parts twice, the real part of the optical conductivity of the many-polaron 
system can be written with a force-force correlation function: 



Re[(T(cij) 



Re / e^-^{[F^{t),FM])dt 
Vhuj'^ mi uo 



(9) 



with F = {i/Ti) [Hq^'P]. The commutator of the Hamiltonian (|l|) with the total momentum 
operator of the charge carriers simplifies to 



N 



(10) 



k i=l 



This result for the force operator clarifies the significance of using the force-force correlation 
function rather than the momentum-momentum correlation function. The operator product 
Fx{t)Fx{'^) is proportional to |Vkp, the charge carrier - phonon interaction strength. This 
will be a distinct advantage for any expansion of the final result in the charge carrier - 
phonon interaction strength, since one power of |VkP is factored out beforehand. Denoting 
Pk = H^=i e^^'^^ , the real part of the optical conductivity takes the form: 



Re[cr(cj)] 



Re I 



k,k' 



Aujt 



Pk'Ok'Vk' + p-k'flk/Vk* 



dt 



(11) 



Up to this point, no approximations other than linear response theory have been made. 



C. LDB canonical transformation for the optical conductivity 



The expectation value appearing in the right hand side of expression (|TT]) for the real 
part of the optical conductivity is calculated now with respect to the LDB many-polaron 
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wave function IV'ldb) @- 

J{k, k') = (V^LDB I [e^''''*/' (pkAkV^k + P-k4^k ) e-*^°*/^ (pk'«k'V^k' + P-i.'a^'V:,)] \ V^ldb) (12) 

(pk'flk'Vk' + p-k'flk'Vk*) f/ 
where U is the many-polaron canonical transformation defined in and H' = U~^HqU 
is the transformed Hamiltonian obtained in ||2|. The canonical transformation of the force 
term is 



(13) 



U ^ (pkCtk^k + P-kflk Vk*) U = pk (Cik - /kP-k) Vk + P-k (flk - /kPk) Vk*. 



(14) 



The terms of lowest order in the electron-phonon interaction amplitude |Vkp are given by 



J(k,k') = |yk|'5kk'(¥'el 

= 2i|Vkp5kk' Im 



e*^'*/Vkake-'^'*/V-k4 - Pkake'^^'^/V-ka^e"^^'*/" 



e^^'*/Vkake-^^'*/V-ka^ 



0) <^ci) (15) 
(16) 



Taking the expectation value with respect to the phonon vacuum, we find 



J(k,k') = 2z|ykr4k' Im{e- 



e^^'*/Vke-^^'*/V-k 



(17) 



This can be substituted in the expression (|Tl]) for the real part of the optical conductivity, 
which becomes 

1 e2 



Re[a{iu)] = -2 



Vhu^ ml 



dt 



(18) 



The right hand side of (|T8D can be written in a more compact form by introducing the 
dynamical structure factor of the electron (or hole) system. 

D. General expression 



To find the formula for the real part of the optical conductivity in its final form, we 
introduce the standard expression for the dynamical structure factor of the system of charge 
carriers interacting through a Coulomb potential. 



-l-oo , 

S{CI,W) = J /v9el 



-oo 



2 ^ 



i^,,)e'^'dt. (19) 



Rewriting expression ( [TB| ) with the dynamical structure factor of the electron (or hole) gas 
results in: 
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Re[a(a;)] = " fc^|y^|25(k, a; - cu^o), (20) 
tiLO rriij J, 

where n = N/V is the density of charge carriers. As noted before, Vk is the electron-phonon 
interaction amplitude and kx is the x-component of the wave vector. Formula (|20|) for the 
optical absorption of the many-polaron system has an intuitively appealing form. 

In the theory of one-polaron optical absorption for weak coupling constants a, the optical 
absorption coefficient as obtained from Fermi's golden rule is ||19|| : 

1 polaron : Re[a{uj)] oc u;^^ ^ A;^! Vk|^(5[^fcV(2mb) - (u; - ulo)]- (21) 

k 

At low densities, the dynamical structure factor S{q, v) is strongly peaked around = u 
and is close to zero everywhere else [^. Substituting a delta-function 5{q^ /2 — v) for the 
dynamical structure factor in formula ( pOD it is easily seen that the one-polaron limit ( pT]) 
is retrieved. The one-polaron result (^Tj) is derived ||TP[ by considering a process in 



which the initial state consists of a photon of energy huo and a polaron in its ground state, 
and the final state consists of an emitted LO phonon with energy TiuJi^o and the polaron, 
scattered into a state with momentum k and kinetic energy (/ifc)^/(2mb) = Ti{uj — cjlo)- 
The many-polaron result, formula (pOf ), is a generalization of this one-polaron picture. The 
contribution which corresponds to the scattering of a polaron into the momentum state k 
and energy Tiioj — culo) is now weighed by the dynamical structure factor S{k,uj — culo) of 
the electron (or hole) gas. 

Formula (|20| ) is reminiscent of the Hopfield formula describing the effect of impurities 



on the optical absorption of metals, which in turn is related to the expression obtained by 
Ron and Tzoar for the optical absorption in a quantum plasma. Formula ( PU] ) also 
represents a generalization of the results obtained by Gurevich, Lang and Firsov These 



authors focused their attention on the many-body effects related to the Fermi exclusion 
statistics, whereas the present analysis will extend the results of to study the influence 
of plasmons and further many-body effects in the system of the constituent electrons or 
holes, as discussed in the next section. 

The advantage of the LDB canonical transformation method 0] for the evaluation of the 
ground state energy of a polaron gas is that the many-body effects are contained in the static 
structure factor of the electron (or hole) system, appearing in the analytical expression for 
the energy. The corresponding advantage of the canonical transformation in the present 
case, for the optical conductivity, is that many-body effects are again incorporated through 
a structure factor, now the dynamical structure factor of the electron (or hole) system. The 
level of approximation made in treating the many-body nature of the polaron system is 
determined by the choice of the dynamical structure factor of the electron or hole system. 

In the present treatment of the electron-phonon interactions, the terms to leading order 
in are automatically taken into account through the variational formulation based on 
LDB 0. As noted before, the use of the force-force correlation function allows to express 
the real part of the optical conductivity as Re[(T(co')] = aJ^{u!,a), so that to lowest order 
in a, Iie[a{uj)] = aJ^{uj, a = 0) where the electron-phonon interaction is no longer present 
in the factor which includes the many-body effects of the electron (or hole) system. A 
possible way to take into account higher-order terms in the electron-phonon interactions 
would be to include electron-phonon coupling effects at the level of the dynamical structure 
factor appearing in JF, in a manner similar to Mahan's treatment of the polaron spectral 



function |24|, or to include multiple-phonon final states in the calculation p5[. 



E. Scaling relation for the optical absorption in two and three dimensions 

The modulus squared of the Frohlich electron-phonon interaction amplitude is given by 



8 



IT/ |2 _ j V V 2mb^LO /99\ 

— ; ^-L^) 

k Ay 2mf,t^LO 

where a is the (dimensionless) Frohhch couphng constant determining the couphng strength 
between the charge carriers and the longitudinal optical phonons, and A is the surface of the 
2D system |2^. In what follows, we will use polaron units {h = mf, = culo = !)• The sum 
over wave vectors in (|20|) can be written as an integral, so that for the three dimensional 
case (with dynamical structure factor Ssd) we find: 

oo 

2 1 f 

Re[cr3D(u;)] = ne^ 3"^^ J ?^'S'3D(g, - u;lo), (23) 



and for the two-dimensional case: 

00 

Re[(T2D(c^)] = ne^ 2'^2^ J (f'S2T>{q,uj - u-lq). (24) 



From these expressions, it is clear that the scaling relation 

Re[a2T>{oJ, a)] = Re[a3Y){i^, 37rQ;/4)] (25) 



which holds for the one-polaron case introduced in ref. [26|, is also valid for the many-polaron 



case if the corresponding 2D or 3D dynamical structure factor is used. 



III. RESULTS AND DISCUSSION 



A. General results 



The expressions (E1,E1) allow us to derive results both for a three-dimensional and for a 



two dimensional polaron gas at T = 0. The choice of a dynamical structure factor for the 
electron (or hole) system allows furthermore to study the different levels of approximation 
(Hartree-Fock, RPA,...) in the treatment of the many-electron or many-hole system. The 
results presented in this section were obtained using the material parameters of GaAs (for the 
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two-dimensional case) and ZnO (for the three-dimensional case). These material parameters 



27|-E9i are summarized in Table I. 



Fig. 1 shows the Hartree-Fock and the RPA result for the 2D many-polaron optical 
absorption spectrum (for GaAs, at a density n = 10^^ cm~^). For reference, the dashed curve 
represents the familiar one-polaron result. In a first step, we discuss the result obtained by 
using the Hartree-Fock expression for the dynamical structure factor of the electron (or hole) 
system in the expressions (|23| , ^^ . 

The Fermi statistics cause the polarons to fill up a Fermi sphere up to kp = [n/(2Tc)Y^'^. 
The optical absorption of the polaron gas resulting from this system is represented by the 
full curve labeled 'Hartree-Fock' in Fig. 1. The spectral weight at frequencies between uji^q 
and 1.4 cjlq in Fig- 1 is reduced as compared to the single polaron case, whereas at higher 
frequencies it is enhanced. A kink appears in the spectrum at u = cjlo + Ep/h, as indicated 
by the dotted vertical line in Fig. 1. 

This can be understood as follows. The absorption process is characterized by an initial 
state consisting of a polaron gas filling up the Fermi sphere (up to energy Ep, at T = 0) 
and a photon with given energy hu, and by a final state made up of an emitted LO phonon 
with energy hui^o and a polaron gas such that one polaron state inside the Fermi sphere is 
not occupied and one polaron state with energy E > Ep is occupied. The incident photon 
can only excite polarons out of the Fermi sea for which huj > huji^o + Ep. A straightforward 
calculation in 2D shows that the fraction of polaron states in the Fermi sphere which can 
interact with a photon of energy Tiuo is given by 

for < culo 
for ..LO < ^ < Ep/h + C.LO • (26) 

h/p 

1 ioT oj > Ep /h + cjlo 

For photon frequencies between culo and culo + Ep/h, the number of polarons which cannot 
participate in the optical absorption process due to the Pauli exclusion principle, decreases 
linearly. At a; = wlo + Ep/h, all polarons can participate. This leads to a kink in the 
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function (p6D describing the number of polarons which can interact with the photon of given 
energy huj, as a function of uj. This is also the origin of the kink in the optical absorption. 
This kink in the 2D many-polaron optical absorption spectrum at huj = TiuJi^o + E-p was 



already noted in [p7 



The full curve labeled with 'RPA' in Fig. 1 is obtained by using the random phase 
approximation (RPA) for the dynamical structure factor of the electron (or hole) system. It 
illustrates the combined effects of the Fermi statistics, discussed in the previous paragraph, 
and screening in the electron (or hole) system. In comparison to the Hartree-Fock curve, 
the main effect is an overall reduction of the spectral weight at frequencies uj > cjlo- There 
is however a second effect, which is the appearance of a contribution related to plasmons - 
this is the subject of the next subsection. 

B. Plasmon-phonon contribution 

The RPA dynamical structure factor for the electron (or hole) system can be separated 
in two parts, one related to continuum excitations of the electrons (or holes) 6*00111 , and one 
related to the undamped plasmon branch |30| : 



SRPA{q,uj) = Api{q)6[uj - ujpi{q)] + Scont{q,^^), 

where ujp\{q) is the wave number dependent plasmon frequency and Api is the strength of 
the undamped plasmon branch [Q. The insets of Fig. 2 depict the regions in the q-v plane 
(g = k/kpiV = m^^uj / {hkp)) where the RPA dynamical structure factor is different from 
zero. The contribution (after substitution of Srpa in (p3|)-(p^)) in the many-polaron optical 
absorption) deriving from the undamped plasmon branch Api{q)6[u! ~ u!p\{q)] will be denoted 
as 'plasmon-phonon' contribution. The physical process related to this contribution is the 
emission of both a phonon and a plasmon in the scattering process. 

Figure 2 shows the result for the optical absorption of the many-polaron gas for the 2D 
case (GaAs, left panel) and the 3D case (ZnO, right panel). For reference, the dashed curves 
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show the one-polaron result. The full curves show the many-polaron results in the random- 
phase approximation. The shaded gray areas indicate the plasmon-phonon contribution. 

Now examine the 3D case (the right panel of Fig. 2). The frequency of the undamped 
plasmon mode lies between ui and U2 where ui = Upi = ^ Anne'^ / rrih is the frequency of the 
plasmon branch at g = 0, and uj2 is the frequency at which the branch of the undamped 
plasmons enters the Landau damping region (whose edge is given hj u = hq^ / {2mh) + 
fik-pq / {2mh)) . The corresponding plasmon-phonon contribution to the optical absorption 
'starts' at cjlo + and 'ends' at cjlo + ^2- These frequencies are indicated by vertical 
dotted lines in the right panel of Fig. 2. 

In the 2D case, the undamped plasmon branch is acoustic-like; for g — 0, uj^i 0. 
Consequently, the phonon-plasmon peak in this case extends from culo up to cjlo + ^2 
where UJ2 is the frequency at which the undamped plasmon branch enters the region of the 
continuum excitations of the 2D (RPA) electron gas. 

In Fig. 3, the evolution of the many-polaron optical absorption spectrum is shown as 
the density of electrons (or holes) is increased. Two effects can be observed for increasing 
density: the reduction of the optical absorption above uj > culo and the shift towards higher 
frequencies of the plasmon-phonon contribution, both in 2D (left panel) and in 3D (right 
panel). The f-sum rule is nevertheless satisfied due to the presence of a central S{uj) peak 
Il3l| in the optical absorption of the polaron gas at T = 0. 



C. Comparison with other theories 



In earlier work, Wu, Peeters and Devreese studied the influence of screening on the 
electron-phonon interaction in a two-dimensional electron gas based on a memory function 
approach using a perturbation expansion in the electron-LO phonon coupling constant. The 



results of this perturbative approach of |2^ for the optical conductivity in 2D, also in the 
RPA framework, are consistent with the results derived from the present method based on 
the variational LDB unitary transformation. These authors found an enhancement of the 
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optical absorption at the frequency where the undamped plasmon branch reaches the region 
of continuum excitations of the electron gas. The present, variational, method extends these 
results by taking into account the entire undamped plasmon branch. 

Recently, Cataudella, De Filippis and ladonisi investigated the optical properties 
of the many-polaron gas by calculating the correction due to electron-phonon interactions 
to the RPA dielectric function of the electron gas, starting from the Feynman polaron 
model and ref. An aspect of the present method is that it is not restricted to the 
random-phase approximation for the treatment of the many-body effects between the charge 
carriers. Cataudella et al. also find a suppression of the optical absorption with increasing 
density. To our knowledge, the plasmon-phonon contribution was not revealed by the work 
of Cataudella et al. [^. 



D. Comparison to the infrared spectrum of Nd2_xCe2;Cu02 

Calvani and collaborators have performed doping-dependent measurements of the in- 
frared absorption spectra of the high-Tc material Nd2-a;Cea;Cu02-y (NCCO). The region of 
the spectrum examined by these authors (50-10000 cm^^) is very rich in absorption features: 
they observe is a "Drude-like" component at the lowest frequencies, and a set of sharp ab- 
sorption peaks related to phonons and infrared active modes (IRAV, up to about 1000 cm~^) 
possibly associated to small (Holstein) polarons |22|. Three distinct absorption bands can 
be distinguished: the 'd-band' (around 1000 cm~^), the Mid-Infrared band (MIR, around 
5000 cm~^) and the Charge- Transfer band (CT, around 10'^ cm^^) p|. Of all these features, 
the d-band and, at a higher temperatures, the Drude-like component have (hypothetically) 



been associated with large polaron optical absorption [p|,|7|,p . 

For the lowest levels of Ce doping, the d-band can be most clearly distinguished from 
the other features. The experimental optical absorption spectrum (up to 3000 cm^^) of 
Nd2Cu02-5 (5 < 0.004), obtained by Calvani and co-workers 0, is shown in Fig. 4 (shaded 
area) together with the theoretical curve obtained by the present method (full, bold curve) 
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and, for reference, the one-polaron optical absorption result (dotted curve). At lower fre- 
quencies (600-1000 cm''^) a marked difference between the single polaron optical absorption 
and the many-polaron result is manifest. The experimental d-band can be clearly identified, 
rising in intensity at about 600 cm~^, peaking around 1000 cm~^, and then decreasing in 
intensity above that frequency. At a density of n = 1.5 10^^ cm~^, we found a remarkable 
agreement between our theoretical predictions and the experimental curve. The experi- 
mentally determined material parameters used in the present calculation are summarized 
in Table I. A background contribution, taken to be constant over the frequency range of 
the d-band, was substracted in Figure 4. The lack of experimental data on several material 
constants leaves us with three adjustable parameters: the electron-phonon coupling constant 
a, the band mass mb, and the density of charge carriers. These parameters were chosen as 
follows: 

• A set of theoretical optical absorption spectra were generated for different val- 
ues of the band mass and densities (for mb = 0.1,0.2,0.5,0.8,1.0,2.0 and n — 
{0.1, 0.2, 0.5, 1.0, 1.2, 1.5, 2.0} x 10^^ cm'^). 

• For each of those spectra, the couphng constant a was chosen so as to fit the tail 

region of the experimental optical absorption spectrum best, using a least squares 
fitting procedure (the tail region is relatively insensitive to many-polaron effects). 

• The best fitting curve, using again a least squares evaluation of the goodness of fit, 
was selected; we found fair agreement with mb = 0.5 me and n — 1.5 10^^ cm~^ at 
« = 2.1. 

This comparison with experiment could not be performed at higher doping content: 
the frequency region of the d-band usually contains a strong non-uniform contribution of 
other optical absorption features (such as the onset of the MIR, the tail of the Drude 
contribution, and the IRAV and phonon modes). To take these other contributions into 
account, additional adjustable parameters would have to be introduced making a comparison 

14 



less convincing. Fortunately, experimental results are available in the form of the normalized 
first frequency moment of the optical absorption spectrum (after substraction of MIR and 
CT band) Re[croxp(cu)] : 



uj 'Re[acy^p{uj)]dijj 

(^) = (27) 

j Re[(Tcxp(t^)]c?c^ 



where cUmax = 10000 cm~^ f^. Calvani and co-workers determined {u) for NCCO samples 
with a varying cerium doping content. Increasing the cerium doping will inject electrons 
in the copper-oxide planes of the material, and increase the 2D charge carrier density in 
these planes. A comparison of this experimental normalized first frequency moment to the 
theoretical one presents the advantage that fewer parameters need to be adapted: only the 
density and the electron hand mass have to be taken from experiment or (if experimental 
values are lacking) fitted. 

The carrier density can be estimated numerically from the effective carrier concentrations 
in the different samples [^] and from a measurement of the two-dimensional Fermi velocity 
performed for one of the samples |l^. As for the other cuprates, the band mass of the 
electrons in NCCO has not yet been determined experimentally |3^, and remains as an 
adjustable parameter. 

Fig. 5 represents the comparison between the present theory and experiment. The 
squares with error bars show the experimental results for differently doped samples of NCCO, 
reported in ||^. The dashed curve shows the normalized first frequency moment of the 
theoretical optical absorption spectrum, integrated over the entire frequency range (tUmax — »• 
oo). The tail region of the many-polaron optical absorption still carries a significant weight, 
just as it does in the one-polaron optical absorption. It is necessary to include the cutoff 
frequency. The full curve represents the theoretical first frequency moment with a cutoff 
frequency tUmax = 10000 cm~\ which corresponds to the experimental cutoff 0. 

There exists a fair agreement between the theoretical and the experimental values of the 
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normalized first frequency moment for the five samples with lowest density, which have a 
cerium doping content of a; < 0.12. These correspond to the squares to the left of the dotted 
vertical line in Fig. 5. For the four remaining samples {x > 0.12) a discrepancy between 
the theoretically predicted first frequency moment for unpaired polarons and the observed 
first frequency moment appears. It has been observed experimentally that the weight of 
the low-frequency component in these samples (with x > 0.12) is significantly larger that 
the corresponding weight in samples with x < 0.12 P|. This was interpreted in P| as a 
consequence of an insulator-to-metal transition taking place around the cerium-doping level 
of X = 0.12. Therefore, it seems reasonable to assume that above this doping level x a 
change in the nature of the charge carriers takes place. One could hypothesize that, as 



the formation of bipolarons is stabilized with increasing density of the polaron gas [34 



bipolarons start playing a role in the optical absorption spectrum. In a variety of other 
cuprates and manganates, the presence of bipolarons has also been invoked to interpret a 
number of response-related properties [BH . 



IV. CONCLUSIONS 



Starting from the many-polaron canonical transformations and the variational many- 
polaron wave function (LDB) introduced in [0] we have derived a formula for the optical 
absorption coefficient Re[cr(ti;)] of a many-polaron gas. We find that Re[cr(u;)] can be ex- 
pressed in a closed analytical form in terms of the dynamical structure factor S{q,uj) of the 
electron (or hole) system, equation (|2^) in 2D and ( ^31) in 3D. In the present approach, the 
electron-phonon coupling and the electron-electron many-body effects formally decouple in 
the expression for Re[(j(u;)]. Therefore, the many-body effects in the electron (hole) system 
can be taken into account by employing any desired approximation to the dielectric response 
(Hartree-Fock, RPA, etc.) of this electron (hole) system. 

In the present work, the dynamical structure factor S{q,uj) of the electron (or hole) gas 
was considered both in the Hartree-Fock and the RPA approximation. 
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The main effect of the Pauh exclusion principle on the optical absorption of the polaron 
gas turns out to be a shift of the oscillator strength towards higher frequencies. This effect 
can be understood in terms of the available initial and final states in the polaron-photon 
scattering process and naturally invokes the Fermi energy Ep of the electron (hole) gas. 

The main effects in the case of the RPA approximation are an overall reduction of the 
optical absorption at frequencies u > culo and the introduction of a novel absorption feature 
which we identified as a plasmon-phonon peak. This plasmon-phonon peak shifts to higher 
frequencies with increasing density, such that a double peak structure can appear in the 3D 
many-polaron optical absorption spectrum, consistent with the observed bimodal polaronic 
band in cadmium oxide ||36[| . 

As a first application of the method presented here, we chose to investigate the optical ab- 
sorption of the interacting polaron gas in the RPA framework. For Nd2Cu02-5 {6 < 0.004), 
similarities were observed (see Fig. 4) between the line shape of the experimental d-band and 
the many-polaron optical absorption as calculated here. To study the density dependence, 
measurements (performed by Calvani and co-workers for a family of Nd2_2:Ce^Cu04_y 
materials) of the first frequency moment of the optical absorption were compared to the 
results of the present theory. We find a fair agreement for the samples with the lowest 
densities (cerium doping x < 0.12). A softening of the first frequency moment of the optical 
absorption for the samples with higher densities (cerium doping x > 0.12) is consistent with 
a change in the nature of the charge carriers at a doping content x = 0.12 inferred in ref. |§] 
from infrared absorption experiments. 
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TABLE 



Table I: Material parameters used in the various figures. The physical parameters for 
GaAs correspond to those of the GaAs-AlGaAs heterostructure ; the material parameters 



for ZnO are taken from [Q. The physical parameters for the neodymium- cerium cuprate 



are taken from ^ and p9[. "n.a." ("not applicable") means that not enough data are 



available to estimate this material parameter. 



material parameters: 




GaAs 


ZnO 


Ndi.ssCeo.isCuOs 


phonon frequency 


huji^o = 


36.77 meV 


73.27 meV 74 meV 


dielectric constants 


£o = 


12.83 


8.15 


n.a. 






10.9 


4.00 


ca. 3. 


band mass 


TTlb = 


0.0657 me 


0.24 me 


n.a. 


coupling constant 


a = 


0.068 


0.849 


n.a. 


polaron length unit 


C^HO = 


5.616 nm 


2.082 nm 


n.a. 


bohr radius 


CtB = 


8.7797 nm 


0.882 nm 


n.a. 



V. FIGURE CAPTIONS 

Figure 1 : The real part of the optical conductivity (proportional to the optical absorp- 
tion coefficient) of an interacting large-polaron gas is shown as a function of frequency, for a 
two dimensional gas (GaAs) from equation |2^. The material parameters are given in Table 
1. The dashed curve represents the one-polaron result, the full curve labeled 'Hartree-Fock 
structure factor' shows the result using the Hartree-Fock approximation to the dynamical 
structure factor of the electron (hole) system, and the full curve labeled 'RPA structure fac- 
tor' is the result in the Random Phase approximation. The dotted vertical line indicates the 
threshold frequency above which all polarons can be scattered into unoccupied final states 
and participate in the absorption process. The broad gray peak in the RPA curve is the 
plasmon-phonon contribution (see also figure 2). 
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Figure 2 : The real part of the optical conductivity is shown as a function of frequency 
for an interacting large-polaron gas in the 2D case (left panel) and the 3D case (right panel). 
The material parameters used here are given in Table I. The dashed curves represent the 
single polaron spectra. The full curve represents the many-polaron spectrum. In this figure, 
the plasmon-phonon contribution to the optical many-polaron spectrum is shown as a shaded 
area. This contribution arises from a process where a polaron, with the absorption of a 
photon, emits a phonon and a plasmon. The inset shows the regions in the q-v plane where 
the dynamical structure factor S{q, v) of the electron (or hole) system used in the optical 
absorption formulas (^ , p^ differs from zero; the Landau damping region and the undamped 



plasmon branch can be distinguished. 

Figure 3 : The real part of the optical conductivity is shown as a function of frequency 
for different densities of an interacting large-polaron gas, in the 2D case (left panel) and 
the 3D case (right panel). The material parameters used for this figure are given in Table 
I. For increasing density, the optical conductivity is reduced. Another effect in the RPA 
approximation is the presence of a peak related to the undamped plasmon branch (see 
figure 2) which shifts according to the plasma frequency. 

Figure 4 : The infrared absorption of Nd2Cu02_5 (5 < 0.004) is shown as a function 
of frequency, up to 3000 cm~^. The experimental results of Calvani and co-workers [^] is 
represented by the thin black curve and by the shaded area. The so-called 'd-band' rises 
in intensity around 600 cm~^ and increases in intensity up to a maximum around 1000 
cm^^. The dotted curve shows the single polaron result. The full black curve represents the 
theoretical results obtained in the present work for the interacting many-polaron gas with 
n = 1.5 10^'' cm~^, a = 2.1 and mb = 0.5 me. 

Figure 5 : The normalized first frequency moment of the optical absorption spectra is 
shown as a function of the density (expressed through the Fermi wave vector). The squares 
represent the experimental results of Calvani and co-workers in a family of Nd2-a;Cea;Cu04 
materials. The dashed curve shows the results from the theoretical two-dimensional many- 
polaron optical absorption, obtained by integrating all frequencies in the calculation of the 
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first frequency moment. Tlie full curve shows the theoretical results obtained by integrating 
up to a cut-off frequency, which is chosen at 10000 cm~^ and which corresponds to the 
maximum frequency in the experiment 0. The material parameters are listed in the figure, 
and the effect of choosing a different electron band mass is illustrated in the inset. The 
points with x < 0.12, to the left of the vertical dotted line, show agreement with the 
theoretical result from the many-polaron theoretical optical absorption, but it is clear that 
the experimental cutoff frequency has to be taken into account. For the samples with 
X > 0.12, a discrepancy between the theoretically predicted first frequency moment and the 
observed first frequency moment is consistent with a possible insulator-to-metal transition 
at X = 0.12 [|. 
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